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The spectral fluctuation properties of various two- and three-dimensional superconducting bil- 
liard systems are investigated by employing the correlation-hole method. It rests on the sensitivity 
of the spectral Fourier transform to long range correlations and is thus an alternative technique to 
study chaotic dynamics. First, we apply the method to the eigenfrequencies which are extracted 
from the measured resonances. Second, we analyze the unfolded raw spectra, including the shape 
of the resonances. The merit of the method lies in a clear separation of the statistics due to the 
positions and due to the shape of the resonances. However, we show that statistical fluctuations of 
the intensities of the resonances have a strong impact on the observable. Therefore, the visibility 
of the correlation hole is studied as a function of the number of independent statistical variables 
entering into the intensities. The visibility improves if independent spectra are superimposed. 
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I. INTRODUCTION 

In recent years, the experimental study of chaotic dy- 
namics in billiard systems has attracted considerable in- 
terest. Unlike nuclei, atoms, molecules or solid state 
probes, billiard systems can be specifically designed to 
investigate certain aspects of chaotic dynamics. This 
"toy model" feature which the aforementioned systems 
lack, makes billiard systems very useful for studies of 
chaotic dynamics. Electromagnetic billiards simulating 
the corresponding quantum systems were experimentally 
investigated for the first time in Refs. ||l|-^. The use 
of superconducting instead of normally conducting bil- 
liards yields an immense improvement in the quality of 
the measured spectra. First results were presented in 
Ref. i. 

The fluctuation properties of a rich variety of systems 
in nuclear, atomic, molecular and solid state physics have 
been studied experimentally and theoretically. In the 
case of fully developed chaos, they are found to be uni- 
versal and very accurately described by Random Matrix 
Theory |mlQ|. Due to general symmetry constraints, a 
time-reversal invariant system with conserved or broken 
rotational invariance is modeled by the Gaussian Or- 
thogonal (GOE) or the Gaussian Symplectic Ensemble 
(GSE), while the Gaussian Unitary Ensemble (GUE) de- 
scribes time-reversal non-invariant systems ^ . Regular 
systems, on the other hand, show significantly different 
fluctuation properties. Remarkably, they are, unlike the 
chaotic ones, not generic and can differ from system to 
system. Although one often encounters the complete lack 



of any correlations, which is referred to as Poisson reg- 
ularity, the extreme opposite, i.e. the totally correlated 
spectrum of the harmonic oscillator, also falls into the 
regular class. 

In order to study the fluctuation properties on short 
and long scales, one commonly analyzes the nearest 
neighbor spacing distribution and the spectral rigid- 
ity [nJlOfl. This requires to extract the positions of the 
levels from the measured spectrum. Thus, the analysis is 
done on the so-called "stick-spectrum" , i.e. the sequence 
of those eigenenergies or eigenfrequencies which could be 
identified in the raw spectra. However, one often has to 
deal with poorly resolved spectra which prevents some 
levels from being found. This "missing level effect" has 
a considerable impact since it counterfeits correlations 
which do not exist. Fluctuation measures that are less 
sensitive to the missing level effect are therefore highly 
desirable. 

In molecular physics, such a technique was devel- 
oped by Leviandier et al. [ni| for the analysis of the 
long range correlations. The properly smoothed Fourier 
transform of the spectral autocorrelation function maps 
the long range correlations onto short scales in Fourier 
space. As compared to fluctuations of regular systems, 
chaotic dynamics causes a considerable suppression of 
this Fourier transform near the origin, a so-called "cor- 
relation hole" . This has been experimentally observed in 
spectra of the molecules acetylene, methylglyoxal and ni- 
trogendioxyd ||ll|,|l^. Recently, by re-analyzing the Nu- 
clear Data Ensemble, it was shown that nuclear spectra 
exhibit the correlation hole, too [pTp^. 
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As pointed out already, fluctuation measures like 
the nearest neighbor spacing distribution or the spec- 
tral rigidity can only be used for the extracted stick- 
spectrum. The correlation-hole method of Ref. [[ll[ , 
however, is also applicable to the raw spectra since the 
Fourier transform separates the statistics of the positions 
of the eigenenergies from the statistics of the intensities 
and widths. A possibly faulty and incomplete extraction 
of the positions and the widths from the spectra can thus 
be avoided. Consequently, the correlation-hole method 
is, in principle, less sensitive to the missing level effect 
and it is worthwhile to study its applicability to other 
physical systems. 

The theory of the correlation hole for realistic spec- 
tra, i.e. including the line widths, was worked out in 
Ref. 1 14 1 in the framework of a scattering model. In 



the same picture, the correlation hole observed by laser- 
spectroscopy in methylglyoxal |ll| was numerically sim- 
ulated in Ref. p5| | . A summary and a qualitative discus- 
sion of the Fourier transform of statistical spectra can be 
found in Ref. |1^. A conspicuous short version of the 
theory of the correlation hole is presented in Appendix 
A of Ref. ||l^. In Ref. |0, the correlation hole is re- 
lated to its classical analogue, the survival probability. 
The correlation-hole method was first applied to billiard 
spectra by KudroUi et al. ||l9[. However, these authors 
studied only the extracted stick-spectra. 

In the present work, we apply this method to the rich 
variety of billiard spectra that have been measured in 
Darmstadt over the last years. We have two goals: First, 
we want to verify the significance of the correlation hole 
for billiards of quite different geometries by analyzing the 
stick-spectra of the extracted levels. Second, we use the 
correlation-hole technique for the unfolded raw spectra 
and discuss the merits and the problems of such an anal- 
ysis. We present a method of superimposing several raw 
spectra in order to improve upon the visibility of the cor- 
relation hole. This is the first time that such a detailed 
study is performed for billiard systems. The excellent 
resolution of the spectra measured in superconducting 
microwave cavities makes the Darmstadt data the ideal 
object for such an analysis. 

After a short description of the experiment in Sec. ||, 
we present the theoretical concepts in Sec. [II. We per- 



form the analysis of our data in Sec. IV and finish with 
conclusions in Sec. M. 



II. EXPERIMENT 

Due to the equivalence of the stationary Schrodinger 
equation for quantum systems to the corresponding 
Helmholtz equation for electromagnetic resonators in two 
dimensions, it is possible to simulate a quantum billiard 
of a given shape with the help of a sufficiently fiat macro- 
scopic electromagnetic resonator of the same shape [0-p| • 

Since 1991, we have experimentally studied several 



two- and three-dimensional billiard systems using super- 
conducting microwave resonators made of niobium. Fig- 
ure |l] shows some of the investigated cavities as well as 
their dimensions. First experiments using superconduct- 
ing instead of normally conducting resonators were per- 
formed in a dcsymmetrized Bunimovich-stadium-billiard 
and a truncated hyperbola-billiard |j,^,^ using the 
2K-cryostats of the superconducting electron linear ac- 
celerator S-DALINAC ^^ in Darmstadt. Recently, a 
desymmetrized 3D-Sinai-billiard |23| , |24|| was experimen- 
tally investigated in a new and very stable 4K-bath- 
cryostat. Note that the electromagnetic Helmholtz equa- 
tion is vectorial in three dimensions and cannot be re- 
duced to an effective scalar equation. Thus, it is struc- 
turally different from the scalar Schrodinger equation. It 
is of considerable interest that the statistical concepts 
developed in the theory of quantum chaos and random 
matrices are also applicable to arbitrary classical electro- 
magnetic wave phenomena. Conceptually, this is similar 
to the study of spectral fluctuations of elastomechanical 
eigenmodes in aluminum 12^J2^ and quartz |g^ blocks, 
which are also described very well by Random Matrix 
Theory. 
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FIG. 1. The geometries of the superconducting resonators 

which were used in the experiments. Note that all systems are 

desymmetrized to avoid parity-mixing. In comparison with 

the full systems this yields a quarter stadium, an eighth of a 

hyperbola and a 48-th of a 3D-Sinai-bilIiard, respectively. In 

the case of the latter, the sides of the drawn cube are 170 mm 

long and the radius of the removed sphere measures half of 

this length, i.e. 85 mm. 

All resonators mentioned above were excited in the fre- 
quency range between and 20 GHz using capacitively 
coupling dipole antennas sitting in small holes on the 
niobium surface. Using one antenna for the excitation of 
the resonator and either another or the same one for the 
detection of the microwave signal, we were able to mea- 
sure the transmission or the reflection spectrum of the 
resonator, respectively, by employing a Hewlett Packard 
HP8510B vector network analyzer. As an example. Fig. g 
shows a typical transmission spectrum of the 3D-Sinai- 
billiard in the range between 6.50 and 6.75 GHz. The sig- 
nal is given as the ratio of output power to input power 
on a logarithmic scale. The measured resonances have 
quality factors of up to Q « 10^ and signal-to-noise 
ratios of up to S/N « 70 dB which makes it easy to 
separate the resonances from each other and from the 



background. As a consequence, all the important char- 
acteristics like eigenfrequencies and widths can be ex- 
tracted with a very high accuracy |gl|,g^,|9[. A detailed 
analysis of the original spectra yields a total number of 
approximately 1000 resonances for the 2D-billiards, hy- 
perbola and stadium, and nearly 1900 resonances for the 
3D-Sinai-billiard. These eigenvalue sequences form the 
basis of the present test of the correlation-hole method. 
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FIG. 2. Transmission spectrum of the 3D-Sinai-billiard 
in the range between 6.50 and 6.75 GHz. The signal is given 
as the ratio of output power to input power on a logarithmic 
scale. The dashed vertical lines mark the extracted eigenfre- 
quencies /^ as used in the stick-spectrum. 



III. THEORETICAL CONCEPTS OF THE 
CORRELATION-HOLE METHOD 



Hence, we study a spectrum of N levels x^^, ^ — 
1, . . . ,N on this unfolded scale. In the following, the level 
number N is always assumed to be large. To simplify the 
theoretical description, we make the further assumption 
that the line shape L{x) does not depend on the intensi- 
ties 2/^ for a given resonance fi in the expression 



I(x) 



N 

E 



y^, L{x - Xf,) , 



(1) 



where the line shape is normalized to unity. The observ- 
able of interest is the decay function, i.e. the modulus 
squared of the Fourier transform of the spectrum 



\C{t)\' - 
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dxl{x) exp(27ria:t) 



(2) 



The Fourier coordinate t defines an unfolded time. This 
expression can be rewritten as the Fourier transform 
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\c{tW 



duj A{uj) exp(27riwt) 



(3) 



of the autocorrelation function 

4-00 

A{uj) = dni{n-uj/2)i{n + uj/2) . (4) 



We summarize, for the convenience of the reader, the 
and some earlier results on the 



basic ideas in Sec. [II A 



two-level form factors in Sec. IIIB. We discuss the the- 
oretical results on the correlation hole for stick- and raw 



spectra in Sees. [II C and HID, respectively. 



A. Basic Ideas 

Consider a spectrum I{E) or /(/) measured as a func- 
tion of energy E or, in our case, frequency /. This 
spectrum is described as a finite superposition of iso- 
lated or interfering resonances of a given universal shape 
L(f) with statistically distributed positions, intensities 
and widths. Although the shape is often known to be 
a Lorentzian p8p9| ], we keep the discussion general. In 
order to analyze the true fluctuations, one removes sec- 
ular variations of the level density, i.e. the Weyl- or 
Thomas-Fermi contribution [pO| p2|. To this end, one 
introduces the smooth part N^'=y\f) of the inte- 
grated level density as the new coordinate by setting 
X — A^Wcyi^j-j £qj, q^j. experimental and theoretical dis- 
cussions. Since we wish to study generic fluctuations, we 
shall henceforth assume that this unfolding procedure has 
been performed and write / = lix). This implies that, 
in the variable x, the mean level spacing is unity every- 
where. 



At this point a comment on the assumed independence 
of the the line shape L(x) of the intensities y^ is in or- 
der. From scattering theory ||3^ it is known that this 
is the case only if the number A of open decay channels 
is very large. Strictly speaking, this is not true for the 
experiments to be analyzed here. The decay channels are 
the antennas that couple the cavity to the external world 
and, hence, A = 2 — 4, see the discussion below. However, 
the assumption of independence is to some extent justi- 
fied here since it affects only the long time behavior of the 
decay function |C(t)p, whereas the correlation hole is a 
feature found in the short time behavior. Moreover, the 
more general theory does also show that the neglect of 
interference effects as done in Eq. (|l|) is justified. Thus, 
it turns out that our simplification does still yield a sat- 
isfactory description of the correlation hole in all cases 
we studied. 

Replacing the energy average in Eq. (0) by the average 
over an ensemble of resonances, one can make use of two 
results of Random Matrix Theory [||J|,|2^ : (i) the y^ are 
independent of the a;^ and (ii) for fJ. y^ v the intensities 
y^^yu are independent of each other. The function (0) 
then can be cast into the form 



A{uj)=Ny^ / dx L{n-uj/2-x)L{VL + ujl2-x) 
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Ny^ / dx' / dx" L{n~uj/2~x') 
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xL{n + u;/2-x")Y2ix" -x') . (5) 

Here, l2(x) denotes the Dyson-Mehta two-level cluster 
function [p|. It describes the two-point correlations, more 
precisely, (1 — Y2{x))dx is the probability to find two reso- 
nances separated by the distance x on the unfolded scale. 
Thus, for large arguments x, the correlations have to dis- 
appear and Y2{x) approaches zero. Note that y^, the 
square of the first moment of the intensities, and j/^, the 
second moment appear in Eq. (|5|). 

We introduce the Fourier transform of the two-level 
cluster function. 



-t-oo 



^2(i) = / dxY2{x) exp{2'Kixt) , 



(6) 



The case that the positions a;^ , Xi, of any two dif- 
ferent resonances fj, ^ v are completely uncorrelated 
is referred to as Poisson regularity |g,0|. Obviously, 
the two-level cluster function must be zero everywhere, 
yPoisson^^;) = 0, and thcrcforc also 6^°'''''°" (i) = which 
results in |C(i)p « Ny'^ — const, even for small values 
of t. There is no correlation hole. 



In the case of fully developed chaos, the statistics of 
the positions is described by the Gaussian Ensembles Q . 
The general symmetry constraints imply that a time- 
reversal invariant system with conserved or broken rota- 
tional invariance is modeled by the Gaussian Orthogonal 
or the Symplectic Ensemble, while the Gaussian Unitary 
Ensemble describes time-reversal non-invariant systems. 
We summarize the results for the form factors Q. The 
situation most commonly encountered is described by the 
GOE, yielding 



and analogously the Fourier transform L{t) of the line 
shape L{x). The function 62 (t) is referred to as the two- 
level form factor Q . Just like the two-level cluster func- 
tion, it has to vanish for large arguments t. Since the last 
term on the right hand side of Eq. dq) contains a convo- 
lution, we can make use of the convolution theorem to 
evaluate the Fourier transform of A{u!). Collecting ev- 
erything, one finds |l^,|lj,|l^,|r^ for non-negative times 



\C{t)\' = Ny^5{t) + Ny^ |L(t)|^ (1 - afosW) 



(7) 



where we denote the ratio of the statistical moments by 



y^ /y^ 



(8) 



The correlation hole is described by the function (l — 

ab2{t)), while |i(t)P describes the decay of the reso- 
nances. Thus, the statistics of the positions is separated 
from the statistics of the intensities. The J-function in 
Eq. (^ occurs since we assumed TV to be very large. For 
a finite number of levels, this contribution will acquire 
a width as discussed in detail in Refs. M,17|. Due to 
the high number of levels in our data, we disregard this 
contribution in the sequel. 

If the resonances are well isolated, then |i(i)P varies 
much more slowly than b2(t). In the extreme case of 
vanishing widths, i.e. L{x) = S{x) the function (M) ap- 
proaches Ny^ for large times t. For the more realistic 
Lorentzian line shape |2g,|2£| one has L{t) — exp(— TrFi), 
where F is the total width, implying that the function (H) 
decays exponentially. 



B. Correlation Hole and Two— Level Form Factors 

For the convenience of the reader, we collect here the 
well known results for the form factors &2(i) introduced 
in Eq. (1). 



lGOE 
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it) = 
1- 
-1 



2t + tln{l 
+ tln{{2t^ 



l)/(2t 



< t < 1 

1)), t>l 



(9) 



This function is displayed in Fig. 0. It corresponds to a 
decay function ([7|) for a — 1 which has been normalized 
by Ny'^. The GSE form factor is given by 



fcGSE 



W = 



t/4+(tln|l-i|)/4, 



0< t< 2 

t > 2 



(10) 



and thus exhibits a divergence at t = 1. Importantly, 
for all fully chaotic cases, we have for vanishing times 
62(0) — 1. Thus, according to Eq. (|^), the deepest point 
of the correlation hole is reached for small values of the 
time t and has the value |C(i)p ~ Ny'^{l — a). 
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FIG. 3. The phenomenon of the correlation hole due to 
the two-level form factor is clearly visible for small values of 
t. The function 1 — &2(i) is displayed for the cases of Poisson 
and GOE statistics according to Eq. (bl). 



C. stick— Spectra 

Suppose that all resonances have a vanishing width, 
i.e. we consider a sequence of levels at positions Xfj. rep- 
resented by i5-functions. If, moreover, all of them have 
the same intensity, y^ = y, Eq. (|]) reduces to 



N 



(11) 



which is nothing but the spectral function. It serves 
as the mathematical definition of the term "stick- 
spectrum" pa. Formally, the probability distribution 
of the intensity y^ is given by a ^-distribution which im- 
plies y^ — y^, hence, according to the definition (ra), we 
find a = 1. Consequently, the function |C(t)P of Eq. (|^) 
becomes 



\C{t)\'^Ny'{S{t) + l-b,{t)) 



(12) 



This is the case of maximum visibility of the correlation 
hole. For a realistic distribution p{y) of the intensities, 
we will always find values of a which are smaller, often 
considerably smaller, than unity. 



D. Correlation Hole for Ra^v Spectra 

The typical raw spectra of our resonators can be de- 
scribed as a superposition of isolated Lorentzian reso- 
nances [ p8|pS| ] with statistically distributed locations, in- 
tensities and widths. This closely parallels the situation 
of isolated resonances in a compound nucleus scattering 
experiment [p3| . Due to a minimized surface resistance in 
the superconducting billiards, the total width is basically 
composed of a sum over partial widths which describe the 
power dissipation into the A decay channels realized by 
the antennas. Thus, each total width is given by 



1 /i — 2_^ 1 Mc 



(13) 



for every resonance labeled by fi and channels denoted by 
c. The small number of decay channels A in our experi- 
ments, i.e A = 2 — 4, leads to a set of strongly fluctuating 
widths r^t . We remark that this allows to examine quan- 
tum phenomena like the non-exponential decay of the 
spectral autocorrelation function [ pl[ . 

The s impli fied theoretical model which was introduced 
in Sec. Ill A, however, is by construction not suited to 
describe this general case but rather explains spectra of 
systems with a constant total width. This situation can 
arise due to a high number of decay channels [^ or due 
to certain peculiarities of the physical system or the mea- 
surement [|ll|,|l^. As mentioned above, the present sim- 
plification concerns the long time behavior of |C(i)p and 
not the correlation hole which shows up on short time 



scales. This simplification is not valid if there are size- 
able correlations between the intensities y^ and the total 
widths F^. Therefore, in the present discussion, it is suf- 
ficient to assume a realistic distribution for the intensities 
2/p. This defines an adequate model for raw spectra and 
has to be viewed in contrast to the (5-distribution of the 



intensities y^ considered in Sec. IIIC. It should be em- 
phasized that we, for reasons of consistency, neglect inter- 
ference terms between the individual resonances. Hence, 
in the terminology of scattering theory, our model applies 
to the case of isolated resonances. 

If the spectrum I{x) has been measured in reflection 
in, say, channel a, the intensities y^ of the /i-th resonance 
are given by 



Vt^ 



r„ 



(14a) 



where F^q is the partial width of the resonance with re- 
spect to the channel, i.e. the antenna, a. Experiments as 
well as the theory of random matrices show that the par- 
tial widths are distributed according to a Porter-Thomas 
law pJfl,Bl|, which describes a fully chaotic system, re- 
sulting in 



a = 1/3 



(14b) 



If, however, the spectrum I{x) has been measured in 
transmission from channel a to, say, channel b, the in- 
tensities y^ are given by the products 



y^ 



J- fia 1 fib 



(15a) 



of the partial widths with respect to the entrance and 
exit channels. Taking T^a and F^;, to be statistically 
independent variables with Porter-Thomas distributions 
one arrives at 



a = 1/9 



(15b) 



These two values of a simply reflect a Gaussian distribu- 
tion for the decay amplitudes which is at the core of the 
Porter-Thomas law. 

Hence, as indicated already, the statistical fluctuations 
of the weight y more or less suppress the correlation hole. 
The (5-distribution considered in Sec. Ill C is much more 



favorable for the visibility of the correlation hole than a 
realistic distribution which will always give a < 1. Re- 
ducing these fluctuations would restore the correlation 
hole. This can be achieved by the superposition of statis- 
tically independent spectra. Suppose that spectra have 
been measured via all, or all possible, combinations of 
the A antennas attached to a given resonator. The posi- 
tions x^ of the resonances are the same in all spectra, the 
intensities y^, however, vary from spectrum to spectrum. 
Thus, the intensities of the spectrum obtained from su- 
perimposing all these spectra will fluctuate much less. In 
the limit of a superposition of inflnitely many spectra, all 
intensities will be the same and we are back to the case 
of the stick-spectrum. 



We want to make this argument more quantitative. 
First, we discuss reflection measurements. Write the /i- 
th intensity of the superposition of all possible reflection 
spectra in the form 



2/m 



/ . -^A'Q 



(16a) 



Under the assumptions that F^a has Porter-Thomas 
statistics and that the average value is independent of 
/i we find 



A 



A- 



(16b) 



For A = 1 one recovers the result ( 14b ) and, as expected, 
this expression approaches unity for large A. 

Second, we turn to transmission measurements. 
Adding up all possible transmission spectra, we obtain 
the ^-th intensity 



2/m 



1 ^ 

9 Z^ r^aT 



fib ' 



(17a) 



,b=l 



Under the same assumptions as above this leads to 
A-1 



A + 7 



(17b) 



Again, this is consistent with the result ( 15b ) for A = 2. 
Moreover, for large A, the expression ( |l7b| ) approaches 
unity, as it should be. 



IV. APPLICATION TO EXPERIMENTAL DATA 



we 



After some general considerations in Sec. IV A 
prese nt th e a nalys is of stick- and raw spectra in 
Sees. IV B and IVC, respectively. 



A. General Considerations 

The Fourier transform C{t) of the measured spectrum 
I{x) does still contain the entire information. As is well 
known, if the experimental data consist of many levels in 
a sufficiently long interval, the Fourier transform can, be- 
fore unfolding, be used to obtain information about the 
periodic orbits of the system |3J| which manifest them- 
selves in a rich structure consisting of many peaks. Here, 
however, we aim at an understanding of the generic sta- 
tistical features of our experimental data. In other words: 
since we are not interested in resolving individual prop- 
erties like periodic orbits, we have to average over all re- 
alizations of the physical system in question. This, how- 
ever, is precisely what Random Matrix Theory did for us 
when we went from the autocorrelation function M) to 



its ensemble average. According to Delon et al. [Q the 
ensemble average can be simulated by applying a smooth- 
ing procedure to the experimental decay function |C(i)p. 
It turns out ||l2| that the most appropriate procedure is 
a convolution of |C(t)p with a Gaussian. Hence, we have 
to compare the theoretical results to the function 



+ 00 



(|C(i)P)= / dt'\Cit')\' 



1 



y/2^: 



exp 



2af 



(18) 



where the variance at was chosen to depend on the time 
as (Tt = t/10. This procedure is referred to as "full Gaus- 
sian smoothing" . 



B. Correlation Hole for Stick— Spectra 

We now analyze the sequences of levels which were ex- 
tracted from the measured spectra, i.e. the stick-spectra. 

In some cases, we perform an additional test of GOE 
characteristics following the discussion of Ref. ||l^. A 
theorem due to Dyson and Mehta pa] states: If one takes 
a GOE spectrum {/i, /2, /a, •■■} and divides it into two 
sequences of odd and even indices, i.e. into {/i, /a, /s, ...} 
and 1/2, /4, /e, ...}, then each of these two spectra obeys, 
after proper unfolding, the GSE statistics. Following this 
idea, the experimentally found stick-spectrum was di- 
vided into two equivalent sets with half the number of 
eigenfrequencies, put together sequentially and the whole 
analysis was repeated. The theoretical prediction is given 
by using Eq. ( p^ ) in the expression (|l^). In contrast to 
the GOE case, the GSE form factor 62(0 has a singu- 
larity at t = 1 due to the different oscillatory structure 
of the two-level correlation function 12(2;). Thus, one 
expects a characteristic peak in the function (|C(t)p) at 
this time t which gives information on correlations on the 
scale of about two mean level spacings. As pointed out in 
Ref. ||l^, the decay function constructed in this way by 
omitting every other level contains information on higher 
than two-level correlations of the original spectrum. 

In Sees. IV Bl and IVB2, we discuss the hyperb ola 
and the stadium-billiard, respectively. In Sec. [VB3, we 



analyze the stick-spectrum of the 3D-Sinai-billiard. 



1. Hyperbola-Billiard 

The upper part of Fig. Q shows the result for the mea- 
sured spectrum of the hyperbola-billiard whose shape is 
displayed in Fig. |l|. The full line is the experimental re- 
sult for (|C(i)|2) according to Eq. (|l|). 

Since the hyperbola-billiard has been proven |]36|,p7| to 
be fully chaotic in the classical limit, one expects pq ] to 
find GOE fluctuations in the spectrum. Indeed, as Fig. ^ 
shows, the agreement with the theoretical prediction of 
Eq. (|l2) with 62 (t) given by Eq. (Bh is very good. Note 



that (|C(t)p) as well as Eq. (12) have been divided by 
NJj^ to allow a comparison with Fig. H. 
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FIG. 4. The function (\C{t)\^) for the stick-spectrum of 
the hyperbola-billiard. Note that the ordinate is divided by 
Ny^ in order to allow immediate comparison with Fig. 3 by 
using Eq. (|l2|). The full line is the experimental and the 
dashed line the theoretical result. 

The test of GSE statistics according to the Dyson- 
Mehta observation mentioned above yields the curve 
shown in the lower part of Fig. 0. Again, the agree- 
ment with the theory is good. Here, to be consistent, 
we have smoothed the theoretical curve, too, by help of 
Eq. (|l8|). We do not exclude the possibility that the dis- 
crepancy between this curve and the experimental result 
ai t — \ might hint at deviations of the higher order 
and the medium range correlations in the spectrum from 
the GOE prediction. Due to the limited amount of data, 
however, we cannot perform more detailed tests which 
would be necessary to make a definite statement. 



2. Stadium-Billiard 

The Bunimovich-stadium-billiard | |39| displayed in 
Fig. is also totally chaotic in the classical sense but 
it has an additional feature: In contrast to the hyper- 
bola, the stadium possesses one neutrally stable and non- 
isolated periodic orbit, the so-called bouncing ball orbit 
which propagates between the two straight and paral- 
lel segments of the geometry. This lends a certain type 
of regular characteristics to the system. Thus, in calcu- 
lating (|C(i)p), we test the influence of this remaining 
regularity on the correlation hole. The result is given in 
Fig. o, where in the lower part the bouncing ball orbit 
has been removed by extracting the term that the bounc- 
ing ball orbit contributes to the smooth part of the level 
density pfl] together with the common Weyl contribu- 



tion in the unfolding procedure. As can be seen from the 
figure, this extraction amounts to a very slight correction 
towards the pure GOE characteristics. Due to the map- 
ping of long range spectral properties onto short times, 
the correlation hole is quite insensitive to the bouncing 
ball orbit. 
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FIG. 5. The function (|C(i)p) for the stick-spectrum of 
the stadium-billiard, as in Fig. 4. The experimental curve 
(solid line) in the upper part includes the bouncing ball con- 
tribution. In the lower part, the bouncing ball contribution 
has been removed. The theoretical prediction (dashed line) is 
the same in both cases. 



3. 3D-Sinai-Billiard 



Finally, we have considered a spectroscopic system that 
generalizes the statistical concepts of quantum chaos in 
the sense that, although it does not represent or simu- 
late a quantum system, its spectral features are found 
to coincide with GOE characteristics. As in the case of 
elastomechanical eigenmodes p5|-p7| , this indicates that 
many classical wave phenomena might follow the predic- 
tions of Random Matrix Theory. The system at hand is 
the 3D-Sinai-biUiard ||J2|,|2|,|^ shown in Fig. |I|. Its 
wave dynamics is described in terms of the vectorial 
Hclmholtz equation in three dimensions. The resulting 
function (|C(i)p) is displayed in Fig. g. As in the case of 
the hyperbola-billiard, we also performed the GSE test 
on our data. 

Note that, in this geometry, the eigenvalues of the 
quantum mechanical Schrodinger equation as well as the 
eigenvalues of the vectorial electromagnetic Helmholtz 
equation show, due to the bouncing ball orbits, slight 
deviations from pure GOE behavior |g^ in other statis- 
tics like the number variance Y? or the spectral rigidity 
A3. It is not possible to draw direct conclusions from 



that to the present system since the ray hmit of the vec- 
torial Helmholtz equation is different from the classical 
limit of the Schrodinger equation. 
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FIG. 6. The function {\C{t)f) for the stick-spectrum of 
the 3D-Sinai-billiard, as in Fig. 4. 

However, it could be speculated that the influence of 
the corresponding bouncing ball orbits in the present sys- 
tem is suppressed due to an effective average over all their 
different lengths. In any case, the effect of the bouncing 
ball orbits becomes visible at comparatively large lengths 
in the spectrum. Thus, in our observable, it cannot be 
extracted with statistical significance since the Fourier 
transform maps it onto values of t which are of the or- 
der of the inverse length in the spectrum. Remember 
that the correlation hole is found at values of t which are 
roughly of the order unity. Nevertheless, a suppression 
of the peak at i = 1 in the experimental curve is seen in 
the GSE test in the lower part of Fig. ||. One might be 
tempted to interpret this as due to those effects. Again, 
more detailed tests would require more data and, in this 
particular case, a thorough theoretical discussion of the 
ray limit of the vectorial Helmholtz equation, too. 



C. Raw^ Spectra 

As far as the stick-spectra are concerned, the phe- 
nomenon of the correlation hole is obviously well un- 
derstood. The agreement with theoretical predictions 
for all considered systems is satisfactory. We now ap- 
ply t he met hod to orig inal and idealized ra w spectr a in 
Sees. |IVCI| and |IV C 2[ respectively. In Sec. |IVC3| , the 
influence of certain statistical fluctuations on our findings 
is discussed and demonstrated using synthetic spectra. 



1. Original Raw Spectra 

We analyze original raw spectra of the hyperbola with 
a total number of A = 3 antennas. Besides the unfolding, 
no further preparation has been performed. On Fig. M, 
the function (|C(i)p) is shown for raw spectra of the hy- 
perbola billiard. The ordinate is in principle as in Figs. 3, 
H and ^: The function (|C(i)p) is given in units of Ny^ 
in order to allow immediate comparison with Fig. and 
all other similar figures, see Eq. (0). It is, however, not 
obvious how to obtain Ny"^ from a spectrum unless one 
identifies and analyzes all the resonances. We proceed as 
follows. 

If one assumes the line shape L(x) to be a Lorentzian 
with width F, then the integral over the square of the 
spectrum I{x) of Eq. (|l|) is 



47r ^ 



47r 



"T3" 



(19) 



This quantity is obtained by numerical integration. The 
width F is obtained by a fit-procedure from the exponen- 
tial decay of (|C(i)p) for t > 1 as is illustrat ed by Fig. 0. 
See also the discussion at the end of Sec. Ill A. From 
these two pieces of information, Eq. ( |l9| ) yields Ny'^. 
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FIG. 7. The function (|C(f)p) for original raw spectra of 
the hyperbola-biUiard. The two parts of the figure show the 
resuh from a spectrum measured in reflection (Sa) and one 
measured in transmission (Sij). The scale of the ordinate is 
discussed at Eq. ([19|). 

In the upper part of Fig. 0, the function (|C(i)p) is 
shown for a raw spectrum of the hyperbola measured 
in reflection mode. Unfortunately, this function cannot 
be interpreted in terms of the correlation hole because 
it is dominated, up to t « 2 by a background which is 
typical for the reflection measurements. Since the radio 



frequency cables we used for the transmission between 
source and resonator are interrupted by certain non-ideal 
cable connectors, the measured frequency spectrum is 
modulated by the signal that is generated through reflec- 
tions at these connectors and the resonator; see the oscil- 
lations of the spectral background in Fig. 1 of Rcf. [£lj . 
Furthermore, the size of the resonators and the lengths of 
the cables are of the same order. Therefore, the period of 
this spectral modulation is of the order of the mean level 
spacing and the power spectrum between t = and t « 2 
is dominated by this artificial background peak. Note, 
however, that this does not at all preclude the analy- 
sis of every given resonance, because F is very small as 
compared to the period of the above oscillations. 

The transmission spectra are free from that problem. 
Nevertheless, the correlation hole is not visible in the 
lower part of Fig. 0, where (|C(i)p) is given for a raw 
spectrum of the hyperbola measurement in transmission 
mode. At the place of the correlation hole, one ob serves 
fluctuations that will be discussed in Sec. IV C 3 below. 



They are due to the statistical fluctuations of the inten- 
sities y^ Note that the value of a = 1/9, expected from 
Eq. ( 15b ) for this case, leaves littl e hop e to see the cor- 
relation hole. According to Sec. [IIF)| , the value of a 
should improve if several spectra are superimposed. We 
have done so for all transmissio n sp ectra (A = 3) of the 
hyperbola. According to Eq. ( |l7b| ), one then expects 
a = 1/5. Figure || shows that this is not enough. As 
compared to the lower part of Fig. H, the fluctuations are 
more suppressed, but the hole itself cannot be identified. 
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FIG. 8. The function (|C(t)p) for three superposed trans- 
mission spectra of the hyperbola. The scale of the ordinate is 
discussed at Eq. (p_9|). 



2. Idealized Raw Spectra 

In order to acquire a deeper understanding of the sta- 
tistical effects that are important for the correlation-hole 
method, we study "idealized raw spectra" . They are 
obtained by providing the experimentally found stick- 
spectra with realistic intensities: We use the experi- 
mentally determined spectral properties of the stadium- 
billiard with A = 3 channels, i.e. the full set of parame- 
ters for the first 950 resonances including their positions 
Xfj, as well as three sets of partial widths F^c which were 



shown to obey the Porter-Thomas distribution [£1| . Note 
that these partial widths have been separately normal- 
ized to unity for each channel, i.e. (F^c) = 1- The re- 
sulting functions (|C(i)p) for different cases are given in 
Fig. 0. In the upper half of the figure, transmission spec- 
tra Sij for i j^ j are analyzed. In the upper left part, a 
single transmission spectrum is used. This means A = 2 
in Eq. (^TJ) and the theory predicts a = 1/9. In the 
upper right part, three transmission spectra have been 
superposed, i.e. A = 3 and the prediction is a = 1/5. In 
the lower half of the figure, reflection spectra Sa are an- 
alyzed. A single reflection spectrum is used in the lower 
left part. This means A = 1 in Eq. (|l^) yielding a = 1/3. 
In the lower right part, A = 3 reflection spectra have been 
superposed and Eq. ( |l6| ) predicts a = 3/5. Note in addi- 
tion, that (|C(i)p) as well as the function in Eq. (R) were 
divided by Ny^. Obviously, there is a strong deviation 
between the experimental and the theoretical curve in 
the case of the single spectra in the left column of Fig. 0. 
For the superpositions, this deviation is reduced. 

N = 950 
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FIG. 9. The function (|C(i)p) for an idealized raw spec- 
trum of the stadium-billiard, as in Fig. 4. See the detailed 
explanation in the text. The experimental curves are given 
as solid, the theoretical ones as dashed lines. 



3. Statistical Fluctuations of the Squared Intensities and 
their Impact 

To give a qualitative interpretation of this deviation, 
we study the statistical fluctuations of |C(i)p in the long- 
time-limit i — > cxD, where the power spectrum is free of 
effects due to level clustering. We will show that the de- 
viation just observed can be attributed to, at first sight 
unexpectedly, large statistical fiuctuations of the squared 
intensities in the spectra. From Eqs. (|l|) and (||) with 
L{x) = 6{x) one obtains 



N N 

\Cit)\' = J2yl+ ^ 






y^y^ exp (27rz(x^ - x„)t) . (20) 



In the limit t ~* oo and due to the full Gaussian smooth- 
ing, the strongly fiuctuating second term of this expres- 
sion is suppressed. We define this limit as 



N 
t—^oo ^ — ^ ^ 



(21) 



The aforementioned ensemble average in the limit iV 
oo yields the quantity 



X= lim X = 7Vy2 



(22) 



which is exactly the expression we used for the normaliza- 
tion in the previous analysis of the idealized raw spectra. 
Now, to estimate the statistical flu ctua tions of X, we 
have to calculate its second moment X"^ and the relative 
standard deviation 



8^c\X 



ix^-x 



X 



(23) 



In the case of a single transmission spectrum 5*^^ with 
A = 2, one has y^ = r^aF^fc which yields 



'^'■'='^ - 7^ \ 



£1 



- 1 



(24) 



This behaves, as expected, like l/\/N. However, since 
the partial widths are Porter-Thomas distributed, the 
higher order moments of F^ obey 



r'= = (2fc-i)!!r 



(25) 



implying that the statistical fluctuations strongly in- 
crease with the order of the statistical moment due to 
the factor {2k — 1)!!. In the present case of A^ = 950 
resonances one obtains dre\X « 0.38 which is in good 
agreement with the experimental curve given in Fig. M. 
For an increasing number of channels A, this relative vari- 
ance shrinks, since the superposition reduces the fluctu- 
ations in the intensities y^. For vanishing fluctuations 
in the intensities, one has Si-dX = 0. Since the param- 
eter a is generated through the statistical moments of 
the y^, both, the normalization as well as the correla- 
tion hole itself, approach the theoretical prediction only 
as the number of included open channels is increased. 

Remarkably, a superposition of reflection spectra even 
for the small number of only A = 3 open channels allows, 
in principle, a re-observation of the correlation hole since 
we have a = 3/5. For the superimposed transmission 
spectra with a — 1/5, the hole is still weak. This explains 
why the hole could not be observed in the superposition 
of the original transmission spectra of the hyperbola in 
Fig. ^. Unfortunately, the favorable behavior of the re- 
flection measurements cannot be exploited because of the 



experimental artifact discussed in Sec. IV C 1 



In order to demonstrate the influence of the number 
A of open channels on the correlation hole and on the 
normalization, we have finally calculated several super- 
positions of synthetic spectra for the hyperbola. We used 



the experimental spectrum of the positions x^ and up to 
A = 10 numerically simulated spectra of Porter-Thomas 
distributed partial widths F^ with (F^c) = 1- T hus, in 



contrast to the idealized raw spectra of the Sec. [VC 2 



where the amplitudes were taken from the measurement, 
the synthetic spectra contain numerically generated am- 
plitudes. Figure Eg displays the results for A = 2, 3, 5 and 
10. Since we treated the spectra as measured in transmis- 
sion mode, we have, according to Eq. ([l^), ratios of the 
moments given by a = 1/9, 1/5, 1/3 and 9/17, respec- 
tively. Obviously, the systematic increase in A leads to a 
continuing reduction of the statistical fluctuations of y^, 
described by 6reiX. For A = 10 the agreement between 
the experimental and the theoretical curve is quite good. 
This demonstrates, for transmission and reflection mea- 
surements, how the correlation hole can be observed by 
superimposing independent spectra which, individually, 
do not show a significant correlation hole. 
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FIG. 10. Dependence of the visibility of the correlation 
hole on the number A of open channels for A = 2, 3, 5 and 10. 
The scale of the ordinate is as in Fig. 4. 



V. CONCLUSION 

In the present work on the correlation-hole method, 
we have evaluated the decay function (|C(i)p) of stick- 
spectra as well as that of original and idealized raw spec- 
tra for various two- and three-dimensional billiard sys- 
tems. 



10 



The stick-spectra lead to the correlation hole as ex- 
pected from Random Matrix Theory. The non-generic 
features implied by the presence of bouncing ball orbits 
in certain billiards did not affect the correlation hole in 
any observable way. 

The fluctuations of the resonance intensities that are 
present in the raw spectra have a strong impact on the 
results. They decrease the visibility of the correlation 
hole. At the same time, they introduce fluctuations such 
that — even with approximately 1000 resonances in the 
spectrum — the decay function (|C(t)p) may fall quite 
far from its expected shape. This again precludes the ob- 
servation of the correlation hole. It is, however, restored 
if sufficiently many spectra with statistically independent 
intensities are superimposed. 
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